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We consider Laplacian Growth of self-similar domains in different geometries. Self-similarity 
determines the analytic structure of the Schwarz function of the moving boundary. The knowledge 
of this analytic structure allows us to derive the integral equation for the conformal map. It is shown 
that solutions to the integral equation obey also a second order differential equation which is the one 
dimensional Schroedinger equation with the sinh^^-potential. The solutions, which are expressed 
through the Gauss hypergeometric function, characterize the geometry of self-similar patterns in a 
wedge. We also find the potential for the Coulomb gas representation of the self-similar Laplacian 
growth in a wedge and calculate the corresponding free energy. 
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I. INTRODUCTION 



The Laplacian Growth (LG) phenomena has been studied for a long time in various physical systems ranging from 
solidification to bacterial colony growth jl|. Displacement of viscous fluid by an invicsid one (to be referred as oil 
and air) between two closely spaced horizontal plates (the Hele-Shaw cell) is the simplest process of this kind. The 
Laplacian Growth is defined here as the dynamics of the oil/air interface (Fig. [1^). For a review, see, e.g., 0, [1]. 
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FIG. 1: Laplacian Growth in a Hele-Shaw cell for the radial a), channel b), and wedge c) geometries. 

The process is unstable and produces in a long time an asymptotical shape which is either a finger [J, Q (in channel 
and wedge geometries) or a fractal with the same universal characteristics [1] as in diffusion- limited aggregation 0, [1] 
(in all geometries). 

It has been recently discovered that the Laplacian Growth is 

• described by the 2D Toda intcgrablc hierarchy Q, 

• equivalent to the evolution of the support of eigenvalues in certain ensembles of random matrices 0, , 

• relevant to the 2D fermions in the magnetic field [ll|, and 

• equivalent to Whitham equations of the soliton theory . 

The most well studied geometry is a rectangular Hele-Shaw cell (Fig. [Jd) where air enters the channel from the left 
side. The boundary dynamics in this system was considered first by Saffman and Taylor in 1958 [4]. They observed 
a relatively stable finger (known today as the Saffman- Taylor finger) formed after a competition of many fingers 
generated by an interfacial instability. They also found a continuous family of traveling wave solutions (in the absence 
of surface tension) and posed a problem of selection of a proper member of the family to describe the observable 
finger. (See the development of the selection problem in |13l|-|la). In the radial geometry (Fig. [T^) an asymptotic 
shape is a fractal with nontrivial universal self-similar features [y], which are still not explained. Another geometry 
which has been studied both theoretically and experimentally is the wedge geometry (Fig. [T];). In this setting air 
enters the wedge filled with oil through the wedge corner and pushes oil outward. A number of exact solutions were 
found in this geometry by different methods [5], 16]-[21]. Some of them are analogs of "traveling wave" solutions. 
They are also called fingers. 

The latter experimental setting is of the great interest for several reasons. One of them is that traveling wave 
solutions in a wedge geometry are the only known explicit solutions of LG described by conformal maps (from the 
unit circle), which derivatives contain non-integer singularities. Moreover, they have a very important property of 
being self-similar. 

It is the latter property that we are going to study in detail in this paper. We prove that the only smooth self-similar 
solution in the radial geometry is an ellipse. The family of self-similar solutions with a singularity at the origin is 
much more interesting. We show that the exact solutions in the wedge geometry first obtained in ^16i] belong to this 
class. 

To study self-similar solutions of general type, we develop a new method which allows us to represent the results 
for the wedge geometry in a more explicit and transparent form. We will show that the self-similarity completely 
determines the analytical structure of the Schwarz function of the moving boundary. Using the knowledge of the 
analytical structure we will derive an integral equation for the conformal map from the exterior of the unit circle to 
the domain filled by oil. We will establish then the equivalence of this integral equation to the Schrodinger equation 
with a well-studied (Poeschle- Teller) potential [l^l which was derived for LG in wedge earlier [l^, [l3| . We will also 
calculate the free energy of the 2D Coulomb gas representation of LG for the case of self-similar solutions in a wedge. 

The paper is organized as follows. In Section [TTl we introduce the notion of the Schwarz function of a curve. In 
Section IIIIl we formulate the Laplacian Growth in terms of the Schwarz function of the moving boundary. Section 
IIVI treats the property of self-similarity in terms of the Schwarz-function. Sections |V] and IVII are devoted to the 
self-similar solutions in the wedge. They contain geometric interpretation, derivation of both integral and differential 
equations for these solutions, proof of their equivalence, and the analytic continuation of the solutions. In section fVlIi 
we present the 2D Coulomb gas representation of the self-similar solution and, finally, in section rVIIII we investigate 
different limits and special cases which are of particular interest for physics. 
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II. SCHWARZ FUNCTION 

The Schwarz function (23j plays a very important role in what follows, so we devote some space here to introduce 
this object. Let a curve C in the plane be given by an equation 

F(x,y) = 0, (1) 

where x and y are cartesian coordinates in the plane. Replacing them by complex coordinates z = x + iy, z — x ^ iy, 
we rewrite the equation for the curve in the form 

G(.,.) = ^(i±^,i^)^0 (2) 

Solving this equation with respect to z (always possible locally for analytic curves) we find 

z = S{z) (3) 

The so defined function S{z) is called the Schwarz function of the curve C. The equation ^ is complex and thus 
contains two real equations, which are mutually dependent, since each of them bears the same information as the Eq. 
(1). The constraint on S is clear after taking the conjugate of equation 

z = SiSiz)). (4) 

This property of the Schwarz function is called unitarity . The Schwarz function is holomorphic and has a great 
advantage over a conformal map description of curves, as it does not imply any particular parametrization on the 
curve. 

The geometrical meaning of the Schwarz function becomes clear if one imagines t hat t he curve C is a mirror. Then 
for the object at a point z close enough to the mirror the image is at the point S{z). For instance, the Schwarz 
function for a circle centered at the origin is S{z) — R'^/z, where R is the radius of the circle. 

Given the Schwarz function T{w) of a closed curve T in the ("mathematical") plane w and the conformal map 
f{w) from the outside/inside of that curve onto the outside/inside of a curve C on the ("physical") plane z, one can 
easily find the Schwarz function of the curve C: S{z) = f{T{f^^{z))), where denotes the inverse function, and we 
adopt the standard notation f(z) = f{z). In particular, if the curve T in the mathematical plane is the unit circle, 
we find that the Schwarz function is given by 

S{z) = f{l/w), z = f{w) (5) 



III. LAPLACIAN GROWTH AND THE SCHWARZ FUNCTION 

In this section we formulate the Laplacian Growth in terms of the Schwarz function of the boundary following [2^ . 

For a quasistatic motion of a viscous fluid between two closely spaced parallel surfaces, the Navier-Stokes equation 
is reduced to the Darcy's Law [l^ which states that the 2D velocity v of the fluid is proportional to the gradient of 
pressure p. In properly rescaled variables it has the form 

V = -Vp (6) 
Ap = (7) 

where the last equation ([7]) follows from the incompressibility condition divv = 0. In the air (zero viscosity) domain 
pressure is constant, p = const. Neglecting surface tension, one concludes that p(x, y) is a harmonic function in the 
oil domain constant along the boundary. In this paper we consider only a single sink placed at infinity. In this case 
pressure p behaves at infinity as 

p{x,y) ~ -^\ogR, R-^oo, (8) 

where R^ — x^ + y^ and Q is the pumping rate (that is the decrease of the oil area per unit time). 

Introducing the holomorphic complex potential W{z) = —p + iip and the complex velocity v = Vx + ivy, one can 
rewrite equation © in the form v = dzW. In the case of a single air bubble in the cell the complex potential is given 
by 

W{z,t)^^\ogw{z,t) (9) 
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Here w{z,t) is (time dependent) conformal map from the oil domain onto the exterior of the unit circle normalized 
in such a way that w{z,t) ~ z/r{t) as \z\ — > oo with a real positive coefficient r(<). This coefficient is called the 
(external) conformal radius of the air domain. In what follows we set Q — tt. In these units, area of the growing 
domain is equal to nt. 

Let us describe the oil/air interface by its time-dependent Schwarz function S{z, t). If at time t a point z is at the 
boundary, z — S{z), then at time t + dt the point z + vdt belongs to the new boundary, i.e., z + vdt — S{z + vdt, t + dt). 
Therefore, on the boundary it holds 

v = S + S'v, (10) 

where dot and prime denote partial derivatives with respect to t and z. 

The unit tangent vector to the boundary equals to dz/dl, where / is the arclength dl — V dzdz — dz^fS' , so 

dz/dl ^ (11) 



In the case of zero surface tension pressure is constant along the boundary, and so its gradient is orthogonal to it. 
In the complex notation we have 

Re^ = 0. (12) 

From unitarity ^ we find S' = \/S' . Together with the previous equation this gives 

v + vS' = Q. (13) 
Using equation (jTU]) and ^ (with Q = tt) we then find 

S^2v = 2d,W = logw (14) 



After analytical continuation from the boundary, this equation is valid everywhere in the oil domain. Equation (jl4p 
connects the change on the Schwarz function of the boundary with the solution to the Dirichlet boundary value 
problem for the oil domain with given sources/sinks. 
One can represent the Schwarz function in the form 

S{z) = S.iz) + S+{z) (15) 

where S- is regular in the oil domain, and S-^- is regular in the air domain. The function S^{z) is given by the Cauchy 
integral 

„ , , If z'dz' 

S+{z,t)^ — f (16) 

where D is the air domain and the point z is assumed to be inside it. (A similar Cauchy integral, with the point z 
outside D, can be written for 5_.) As the function dzW has no singularities in the oil domain, we find 

S+iz) = (17) 

From this we conclude that all harmonic moments of the oil domain 

tk^T^i zz-'^dz-^ (f S+{z)z-''dz, k>l, (18) 
are integrals of motion. This fact was observed by S.Richardson [2^. 



IV. SELF-SIMILAR SOLUTIONS. GENERAL RELATIONS 



Self-similarity expresses the fact that the evolution of the domain is just a dilatation of its shape. It means that if 
initially the Schwarz function of the boundary is S{z), then at any other time z/r(t) — S{z/r{t)) for the points on 
the boundary, where the real dilatation parameter, r(i), initially equals to 1. In other words. 



S{z,t)^r{t)S{z/r{t)). 



(19) 



FIG. 2: A singular self-crossing domain. The thick lines are the branch cuts of the function S+. Spacings between two adjacent 
fingers are called fjords. 



Let us take r{t) to be the conformal radius of the growing (air) domain. It is defined as the derivative at infinity of the 
conformal mapping from the exterior of the unit circle onto the oil domain normalized by the conditions that infinity 
is sent to infinity and r{t) > 0. Then S{z) is the Schwarz function of the oil domain scaled so that its conformal 
radius is 1. 

The self-similar reduction of ea. p^ reads 

S{z) - zS'{z) ^2Ad,\ogw{z) (20) 

where r^(i) = t/A and w{z) is the conformal map of the domain with conformal radius 1. Taking the contour integral 
of equation (|20p over the boundary we find that the constant A is the area (divided by tt) of the domain with conformal 
radius 1. 

Since the r.h.s. of (j20p is an analytic function in the exterior of D, this equation implies the relation 

S+{z) - zS'+{z) ^ 

with the general solution 

S+{z)^sz (21) 

where s is a constant. This solution corresponds to an ellipse with the axes ratio (1 + s)/(l — s). We thus see that 
regular self-similar solutions are exhausted by ellipses. 

This result is not surprising and can be seen right from equations (|17[) and (|18[) . Under the dilatations z — > rz the 
fcth harmonic moment changes as tk — > tkr^~'^ . On the other hand, the harmonic moments are integrals of motion 
and as such cannot change. This leaves as with the only possibility to have all the moments equal to except ^2, 
which is exactly the statement (I?!]) . 

However, there is a family of singular self-similar solutions of equation (|20p for which the interface may cross itself 
at the origin, see Fig. [2] (to be prepared for the next section, the domain looks like ZAr-symmetric one but in general 
no additional symmetry is implied). For such solutions the constant s in (PT|) is different in different sectors of the 
plane and, therefore, S+ has jumps across some rays with argz = const starting at and coming to oo. We stress 
that these rays are singularities of S'+, so they do not depend on t. Let F be such a ray. Take the jump of the both 
sides of eq. (PH)) on this ray. Since the r.h.s. is regular, we get 

[S{z)]r-z[S'{z)]r^Q, zeT 
(the jump is defined as [S'(z)]r — S{ze^^°) — S'(ze~*°), z £ F), hence 

[^(z)]r = c(r)z, zer 

The jump of S{z) on the branch cut is related to the geometry of the interface around the origin. Let us choose 
coordinates so that the branch cut is along the positive real semi-axis. Consider a vicinity of the point z = 0. The 
bmction S'{z) has a constant jump on the positive part of the real axis equal to c c{r). On the other hand, we 
know that the unit tangent vector to the boundary curve is pT|) l{z) = 1/ ^JS'{z), so using the notation in Fig. [31 we 
write = = 5"(0+), l^"^ = e-^^C'^-^^) ^ s'{Q-) (here 0± means ee±*° as e 0) and so the jump is 

c = 5'(0+) - 5'(0_) = e~2''9i _ e2»«2 



FIG. 3: The jump of the Schwarz function is related to the fjord angle. 



Introducing 9 = 9i + 62 (the fjord [s^l angle) and = 0i — 6*2, we obtain the desired relation 

c= -2ie""^sin 6* (22) 

We see that if c is purely imaginary, then the fjord is symmetric (0 = 0) with respect to the real axis. 

The connection between the jump of 5+ and the fjord angle can be understood purely geometrically. Let us consider 
a symmetric fjord, i.e., the ray F is the symmetry axis. Very close to the origin the boundaries are almost straight 
lines. Let us consider a point ze+'°, where z is on F. Its reflection in the boundary lies at ze*^ = S{ze^^^). For the 
point ze^'" we have ze^^^ = S'(ze~*°), so 

[Siz)] = 5(ze+'0) - 5(ze-*) = (e-'*' - e'^)z = cz, 

from which equation (|22|) follows. 

Using the representation of the Schwarz function via conformal map ([5]), equation (|20p can be rewritten in the 
familiar form 

- 2A 

f{w-')d^fH - fiw)d^fiw-') = — . (23) 

w 

The constant A can now be expressed as 

A=ir-<f fiw-')df{w). (24) 

We note that equation ((23l) looks as the Wronskian of f{w) and f{w~^). This fact suggests that they are two 
linearly independent solutions of a second order differential equation of the form 

w^f'iw) + wf'{w) + V{w)f{w) = 

with a real "potential" V, such that V{w) = V{w^^). Below we shall see that this is indeed the case. 

V. SYMMETRIC SELF-SIMILAR SOLUTIONS AND SOLUTIONS IN THE WEDGE 

In this section, we study a family of self-similar solutions with the ZAr-symmetry (Fig. [2]). Equivalently, these are 
self-similar solutions in the wedge with the angle — 2tt/N. 

Let io — e^'^*/^ be the primitive root of unity of degree N. The Z^y-symmetry means that if z G D, then ujz G D. For 
the conformal map this mean that f{u!w) = ojf{w) and for the Schwarz function that S{ljJz) = uj^^S{z). Moreover, 
from the integral representation pop we see that the latter property holds for 5*+ and 5_ separately. 

We need to introduce the following notation: Let Sk be the sector ^(fc — 1) < argz < and let Fk — I^k+i <^^k 
be the ray argz = ^k, see Fig. |31 
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FIG. 4; Thick lines Tk are the branch cuts in 5+. Enumeration of sectors Sfc is also shown. 

A. The integral equation 

Combining self-similarity {S^{z) ~ Sf^z if z G Sk) with Zjy-symmetry (S'_|_(ajz) = u!~^S-i-{z)) we find that s^+i = 
Lu~'^Sk and so 

Sk = Lu-^'^so (25) 

Therefore, 

[^(.2)]r, =cc^-''z, zeTk (26) 

where we have set c = (t<^ — l)so- These jumps are singularities of 5'-|-[3l|. 
The definition of the Schwarz function written in the form 

/» = Sifiw-')) (27) 

implies that singularities of f{w) inside the unit circle exactly correspond to singularities of S{z) outside the domain 
D, i.e., to singularities of S+{z). In our case f{w) = f{w) due to the reflection symmetry with respect to the real 
axis, and singularities of /(w) are jumps on N intervals argui = ^/c, < juij < 1 inside the unit circle, which we 
denote by 7^. Therefore, we can write f{w) as the Cauchy type integral 

(the integral is over 7^). From ([?7|) we conclude that 

[fiw)]^, - - [^+(z)]r„_4^^(„_,) = -cco-'^fiw-') (29) 

(The minus sign is due to the interchange of the two sides of the branch cut when one inverts w [33|.) Plugging this 
into (|28p . one obtains, after some transformations: 



N-l i-oo 

fc=o 1 ^ 



Since f{uj''x) — w'^/(a;), the sum is X^fc"^ '^/(l ~ oj'^xw) which is easily calculated to be Nxw/[\ — x^w^), and so 
we obtain the integral equation for the conformal map: 



wf{x) dx 



(30) 



with the "boundary condition" /(I) = 0. Let us denote the fjord angle hy 9 — 7t/3, then c ~ — 2jsin 7r/3 (from eq. 
(|22| ) and the integral equation acquires the form 



TV sin 77/3 wf{v)dv 
f{w)=w + / Y^^^, /(1) = (31) 



Without using eq. (P^ . one can find the coefficient in front of the integral by noting that the fjord angle 7r/3 implies 
that /(w) ~ (ui — 1)^ near the point w = 1. Expanding eq. pop around this point, one obtains c = — 2zsin 7r/3. 
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B. Geometrical interpretation 



There is another way to obtain the Equation (|30[) . namely, from the geometrical arguments applied to multi- 
logarithmic solutions of 2D LG j23| . 

We start from a simple observation that each logarithm in the standard multi-logarithmic solution (anzats) for the 
conformal map of the LG process asymptotically describes a single "U" shaped fjord, a fjord with parallel walls. Let 
us consider such an anzats: 

f{w) = rw + '^aklog{l - ak/w) (32) 

k 

where afe's are constants of motion, while all {\ak\ < 1) and r are functions of time. Asymptotically at large time 
r — > oo and \ak \ 1. Integrals of motion for such a solution are 

Pk = /(l/ofe) = r/cLk + ^ ak' log (1 - ak'CLk) (33) 
fe' 

At large time, when all akS are close to 1, we can neglect all logarithms except one. Imagine for simplicity that ag is 
real and oq — *■ 1, then we can write 

/3o w r/ao + ao log(l - al) Ki r + ao log(l - ao) -I- ao log 2, (34) 
/(I) w r + ao log(l - ao) = - log 2. (35) 

This means that asymptotically the bottom of the fjord is at the point Po ~ «o log 2. 

Let us now take a number 5 such that (5^1 and (5 ^ 1 — ag, then using e'* 1 -f i(5 we can calculate 

/(I + i(5) w r + ao log(i(5), 
/(l + i(5)-/(l-i(5) = maQ. 

This shows that the width of the fjord is constant and equals Trjaol. The direction of the fjord is given by the phase 
of ao. 

Using this simple correspondence between the logarithmic singularities and the "U" shape fjords we can construct 
a so-called "V" shape fjord, i.e., a fjord which walls form a constant angle between them. 

It is clear that for this construction we need to put many identical fjords in such a way that they have equal distance 
between their vertices). All fjord vertices must be on a ray which direction is the same as the direction of the fjord. 
For the Z^r-symmetric case we then write: 

oo N~\ 

/(ly) =ru; + ^ ^ awMog(l-a/wV^«): (36) 

/?, = = /(1/aO, (37) 

where a/, a, A are real and A is the distance between the nearest /?s. The numbering of a/ starts at the largest a, 
which is 1, so 

ao = 1. (38) 

We now want to take the limit of A, a — s- 0. At this limit the spacing between a's becomes small, and we can think 
of a/ as a smooth function a(l\ or alternatively /(a). We then can go from summation to integration and write 



l{a) 



^ /"OO ^ /■ 

/(w) = + a ^ 'J' / log(l — auj^ /w)dl — rw -\- a ^ ljj'^ I 
k=o -^0 k=o 



:da, (39) 



l{a) = ^/(1/a), (40) 



where we integrate by parts and use ([55)1 . Substituting eq. into ([5^ . we get 

" /, 

fc=o 



f{w) = ™ + X H ^ / 



(41) 



from which (|5D|) follows. We note that the relation /?; — lA (see (P7)) ) is crucial for the derivation and it is this relation 
that leads to the linear integral equation. 
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VI. INTEGRAL AND DIFFERENTIAL EQUATIONS FOR THE CONFORMAL MAP 

A. The integral equation. Direct solution 

In the previous section we have derived the integral equation for the conformal map f{w) (see (|3T|) ). It is an 
inhomogeneous equation. Let us transform it into a homogeneous one. Putting w = 1 in (pij). we get 

TV sin ttP f{v)dv _ ^ 



TT 7i -1 
that allows us to rewrite (1511) as 



After the substitution 

f{w) = w{l - w-^fF{w-^) , w-^ = x , K = 2/N 

the equation acquires the form 



(1 - xf-'F{x) ^ — ^ / \ _ J F{t)dt (42) 







which makes sense for arbitrary k < 1. The structure of the r.h.s. suggests to look for solutions in hypergeometric 
functions. Specifically, we use the identity 







(1 ^)-i r(a)r(c-6) f a, c-b, 1 . x 



r{c-b + a) \ c, c-b+a ' x-1^ 

(taken from [2^, section 7.512). Notice that ate— 6 + a= l the function 3F2 reduces to 2-F1 and the r.h.s. becomes 

a, 1 — a X 



X- 1 



Using the identities 



and 



we see that the function 



(l-2;)-^r(a)r(l-a) 2F1 



r(a)r(l-a) = 



solves eq. (j42|) . 

Summing up, we have obtained the explicit formula for the conformal map: 



(43) 



For Z^v-symmetric patterns k — 2/N and 7r/3 is the fjord angle. Equivalently, this solution describes the self-similar 
growth in the wedge with angle = ^tt/N. Being restricted to the sector < argw < ttk, the r.h.s. of ([331) makes 
sense for all < k < 1 and provides the solution for the self-similar growth in the wedge with angle = 
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We note some simple properties of the map (|43l) . First, there is the symmetry /3 ^ 1 — /3. Second, at /3 = the 
hypergeometric function is identically 1 and the function (|43p is simply f{w) = w that corresponds to the self-similar 
growth of a disk or sector. In the opposite limit, the function ((43|) gives the conformal map to the complement of a 
needle-like domain (a degenerate finger). Let us consider two cases separately: i)0<K<l/2 and ii) 1/2 < k < 1 
(the case k > 1 is discussed in Appendix A). In case i), at /? = k the hypergeometric function is again identically 
1 and f{w) = {w^^^ — w^^^'^)'^ is the map to the complement of a needle growing from the origin. In case ii), at 
/3 = 1 — K we use the symmetry /3 — > 1 — /3 and again get the same mapping. Note that in case ii) for values of f3 in 
the range 1 — «</?< 1/2 a critical point of the function (a value of w at which /'(w) = 0) is outside the unit 
circle, so the map (|43|) is no longer conformal. Such values of /3 are thus non-physical. 

To summarize, the physically meaningful range of the parameters is as follows. If k does not exceed i, then /3 can 
be any number between and k, and tt{k — /3) is the angle of the finger (the angle between the tangent lines to the 
finger boundary at the origin). If i < k < 1, then f3 can not exceed 1 — k. As /3 approaches 1 — k, the finger becomes a 
very thin needle. Formally, in this limit the angle of the finger tends to 7r(K — /3) (not to zero!) but the neighborhood 
of the origin where this angle properly describes the form of the finger vanishes as /? tends to 1 — k. 

With the conformal map (|43|) known explicitly, one is able to find some geometric characteristics of the growing 
self-similar fingers. In particular, the distance from the finger tip to the origin is 



(the conformal radius is 1). Using some identities for the hypergeometric function, one can represent it in a more 
explicit form: 

r(^)r(i-«) 

i? = 2''~i ^ " ^ , ^. (44) 

r(/3 - k)t (1 - ^ 



The area of the finger with conformal radius 1 will be found below. The result is ^ttkA, where the constant A (see 
p4)) ) is given by 

^ r(i - .)r(-.) 

r{/3- k)T{1- (3- k)' ^ ' 

For physical values of the parameters k and /3 this constant is always positive. For non-physical values (1/2 < k < 1 
and 1 — K < P < 1/2) formula (|i5)) gives a negative area. 



B. Differential equations for the conformal map 

It turns out that the solution to the integral equation (pij) obeys a differential equation of second order (that is not 
surprising because the solution is expressed through the hypergeometric function). For the first time it was written 
in [13| . Specifically, the equation reads 

wdUwd^f) + / = /. (46) 

It has two linearly independent solutions: 

]{w) = w{\ - w-^'f 2i^i ; w-""^ , (47) 



g{w) = w-\\-w-^f (48) 

The solution f{w) (growing at infinity as w) is our conformal map (j43[) . The meaning of the second solution (tending 
to at infinity) is clarified below. 

It is useful to present the integral and differential equations in another form which may be more suggestive for some 
purposes. Set w = e^^l^ = e*^^ and /(e*^^) = '^[x). In this notation, the equations are: 
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The integral equation: 



V'(^) = e''" + —T^ I i_J. J) - (49) 



" Jo 
• The differential equation: 

i^"{x) + ^fc^ ^(x) = K^i^ix). (50) 
sinh a; 

The latter is the Schroedinger equation in the singular potential V{x) = — We need solutions on the positive 

real axis such that ipi^) = 0. The equation has two linearly independent solutions 'iP±k{x), where 

V-^x) = e«-(l - e-^^f 2F, ^-"^ ; e''^^ . 

The solution that grows at infinity, ip^^ix), corresponds to the confornial map. 

In order to show that ip^ obeys both integral and differential equations, we rewrite them as the eigenvalue equations 
I^p — ^i^T^V'i Dip ~ (3{(3 — l)ip for the operators 

^OO K.{X+X') • U 

ItP(x)= -r- ifj(x')dx', 

Jq sinh [x + x') sinh x' ' 

Dipix) = smh^x{d^ - K^)tP{x). 

A direct calculation shows that the operators /, D commute on the space of functions vanishing at the origin, and 
thus have a family of common cigcnfunctions. The correspondence between the eigenvalues is established by their 
local behavior around the origin. 

A change of the variables leads to another form of the differential equation, where the parameters k and \ — (3 get 
interchanged. Namely, introduce the new variable ^ such that 

cosh^ = cothx (51) 

(equivalent forms: sinh^ ~ — (sinha;)^^, ^ — log tanh(a:;/2)) and set 

^^(0 = (sinhe)^/V(a:(0). 

Then eq. (|50p becomes 

C. Analytic continuation of the conformal map 

An important information is encoded in the analytic continuation of the function f{w) to the interior of the unit 
circle. It can be achieved by means of the identity 

\ c 'V"r(5)r(c-a) ^-'^ a-b+1 

T{c)r{a-b) , fb, 6-c+l 1 

+ r(a)r(c-6) ^~'> '^'[ 6-a + l 

In our case this identity gives: 
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sm-jv \ ^ N 

where 

1 (P — Kji (1 — p — k) 

Taking appropriate branches of the multi-valued functions, we have 

Therefore, we can compactly write the analytically continued function in the form 

f[w)=Bk![w-^) + Ag{w-^), weSk, (54) 

where 

B,^.^^-^'^, (55) 
sm^ 

A is given by ([53|) . and giw) is the second solution (|48|) of the differential equation ((46|) . (Note that f{w^^) and 
(7(w~^) are also solutions of this equation.) 

In fact ((54)) is nothing else than the decomposition of the Schwarz function S = S+ + To see this, we rewrite 
([5i)l in the equivalent form 



which makes it clear ^ that the l.h.s. is S{z) while the first term in the r.h.s. is S+{z) (cf. (HH)). Therefore, the 

second term is just S-{z) = g(f^^{z)). This observation allows us to compute the area of the domain with conformal 
radius 1 and to conclude that it is equal to nA, in agreement with the previously introduced notation ([24|). Indeed, 
we write: 

Area(D) = ^* S{z)dz ^ (f S+{z)dz + ^ (f S-{z)dz . 
2* JdD 2i 2i Jqq 

The first integral vanishes since S'+ is analytic inside D while the second one is 

Y^j 9{w)df{w) ^ ttA 

(by taking the residue at infinity). 

VII. THE 2D COULOMB GAS REPRESENTATION FOR THE SELF-SIMILAR SOLUTIONS 

As it has been shown in [l^, H^, the LG dynamics (with zero surface tension) can be simulated by the 2D 
Coulomb gas in an external field in the limit Af 00, where Af is the number of charges: 

ZAt)= fl[\z:^z,\^l[e^^(^-^')d^z,, (56) 

i<3 1=1 

where the potential W is determined by the "positive" part of the Schwarz function: 

W{z,z) ^ -\z\^ + 2TZe [ S+{z')dz'. (57) 
Jo 
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Specifically, it can be shown that in the limit — > oo the charges uniformly fill a domain D{t) in the plane and the 
shape of this domain depends on t according to the Laplacian growth equation. For the self-similar ZTv-symmetric 
solutions a simple calculation gives 

2 sin ^ V / 

This potential is continuous across the boundary between the sectors and enjoys the ZAr-symmetry. 
The free energy of the Coulomb gas (|56p . normalized in a special way, 

plays an important role. (As is shown in fl3|, it is related to tau-function of certain integrable hierarchies.) To find 
it for the self-similar solutions, we use the integral representation 

2nJ'{t)= [ W{z,z)(fz + t [ log\z\^(fz (60) 

JD(t) JD(t) 



and the relation 



21ogr(t) (61) 



{r{t) is the external conformal radius of the domain D(t)) obtained in [T3|. Since the potential is homogeneous, 
the self-similar ansatz plugged into ([50)1 yields 

2nT{t)=r^{t) [ Wd^z + Ar\t) f logir^it)\z\^)d^z , 

Jd{A) Jd{A) 

where the domain D(A) has conformal radius 1 and r'^{t) = t/A. Therefore, we can write 
and calculate the second t-derivative: 

d^T 2 , t 3 



= _FiA)+log- + - 



df^ A^ ' ' ^ A 2 
The r.h.s. must be equal to 21ogr(i) = log(t/A), hence J^{A) ~ ~jA^ and eventually we obtain 



^(^) = i^2logi-^^^ (62) 



where A is given in ([53 



VIII. SPECIAL AND LIMITING CASES 



The trivial special cases are f3 = (no fjords, the domain is the disk, f(w) = w) and f3 = 2/N (the domain is a 
collection of symmetric needles growing from 0, f{w) = w{l - w;-^)^/^). A non-trivial special case, where the 
solution is available in elementary functions is A*" = 4 (wedge with the right angle). At A^ = 4, the conformal map is 
i 

/(y;) = I [(1 - W-^il + ii;-2)l-/3 + (1 + ,„-2)/5(l _ ^i,-2)l-/3] 

and A= {1 — 2(3) costt/?. A direct substitution A^ = 1 leads to a singularity. To resolve it, we set N = 1 — e, (3 — Aet2 
and tend e ^ 0. The general formulas then give: S'+(z) = 2^2^, f{w) = w + 2t2/w, A = 1 — At\, i.e., the domain is 
an ellipse with the second harmonic moment t2 (in the normalization used in fo']). The case N = 2 gives, in a similar 
way, the same family of ellipses (restricted to the half-plane). See also Appendix A. 
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The most interesting limit is A'' ^ c». In this case we expect the sohitions to reproduce, after a proper rescahng, 
the SafFman- Taylor fingers in the channel geometry. In the rest of this section we show that this is indeed the case. 
First of all, let us change — s- w and consider the conformal map 



in the limit k = 2/A — > and /3 so that (3/k remains a finite parameter. We set 

-A), 

so that ttAk is the angle of the finger (tending to zero in the limit). Second, passing to the channel of a finite width 
(say 2tt) requires the rescaling of the physical plane z Nz. At last, in order to be closer to the finger tip one should 
shift the origin. Taking this into account, it is not difficult to see that the correct limit to the channel geometry is 

/ch(w)- lim [Nh{w)-N] (63) 

Performing this limit with the fi{w) given above, one obtains 

/ch(u') =logu; + 2(l- A)log(l-u.-i) (64) 

which is the standard form of the SafFman- Taylor finger with the relative width A. Note that the hypergeometric 
function behaves as 1 + O(k^) and so does not contribute to fch{w). 

It is also instructive to derive the limiting form of the time-dependent conformal map and to follow in detail how 
self-similarity in the wedge turns into translational invariance in the channel. At nonzero k the time evolution is just 
a similarity transformation, so we can write, taking into account the shift of the origin: 

Uiw,t) = hm f'-^I^NMw) - A^ 

where Tq = AN^, i.e., r{To) = N. Expanding 

r{To) + t / tV^"^ , t 



r{To) V To J 2To 

we get 

tN 

fch{w,t) = /ch(w) + — . 

It remains to rescale time as t — 2NT, where T is kept finite (note that still t <^ Tq), then the last term becomes 
T/A. The limiting value of A is 1 — — A, so we obtain the time-dependent SafFman- Taylor finger in the standard 
form 

fchiw,T) = Uiw,t) = |+logw + 2(l-A)log(l-u;-i). 
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APPENDIX A: SHAPES OF THE DOMAINS 

In this Appendix we discuss the range of parameters of the conformal map (|43p and the corresponding shapes of 
the air domains. It is convenient to adopt the function ([43]) to the case of a single sector in the physical plane: 

f{w) = w^/'{l l/wf 2F1 f - ^ ; 1/u;) . (Al) 
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(we remind that in our notation the wedge angle is 8^ = ttk and the fjord angle is Q fjord = /Stt). 

First, we notice that at /3 = the hypergeometric function is identically 1, so the conformal map is just /(w) = w^^^ 
- the map from the exterior of the unit circle to the exterior of the sector. 

Further investigation should be done separately in the different regions: i) < k < 1/2, ii) 1/2 < k < 1, iii) 
1 < K < 3/2, and iv) 3/2 < k < 2. Two special points k = 1 and k = 2 require separate considerations. 

In the region i) the allowed values of /? are < (3 < k. At (3 — k the hypergeometric function is again identically 1. 
The map then is f{w) = {w^^'^ — w"^/^)" and the air domain is a needle. The conformal map has a critical point (a 
point where f'{w) =0) w = —1 right on the unit circle. If we formally consider the case (3 > k then we see that the 
critical point is outside the unit circle. The map (|A1[) is then no longer conformal. 

In the region ii) the allowed values of j3 are 0</3<l — k. At/3=1 — Kwe use the symmetry /? ^ 1 — /3 to see 
that the map is again given by f{w) — (w^/^ — w^'^/'^Y ^-nd the air domain is again a needle. The critical point is at 
w — —1. If we increase the value of /3, then the critical point moves outside the unit circle and the map is no longer 
conformal. 

In the region iii) the allowed values of (3 are 0</3<k— 1. At (3 = n — 1 the hypergeometric function is just 
1 + l/w. The function (jAl[) is then f{w) — {w^^'^ — w~^^^)'^~^{w^/^ + w^^/'^). This gives a needle-like domain but 
the needle is "two-sided" and rotated by 7r/2 (see the figure). The map has two critical points on the unit circle at 
the points w — exp{zLi(j>o) , where tan^o = ^/k — 1. This corresponds to the two tips of the needle. If wc increase the 
value of /3, then both critical points move outside the unit circle and the map is no longer conformal. 

In the region iv) the allowed values of (3 are < (3 < 2 — k. At (3 = 2 — k we use the symmetry /3 ^ 1 — /3 and 
again see that the situation is identical to the one in region iii). There are two critical points right on the unit circle. 
Both of them go outside should we increase the value of /3 even further. 

All the four cases are illustrated in the figure. 

*:=l/3, 0<j8</t 




The special value k = 1 requires careful consideration. We set k = 1 — S and take the limit S ^ 0, P ^ with 
the arbitrary ratio a = (3/6 < 1. The hypergeometric function is then 1 — a/w and the conformal map (|A1[) becomes 
f{w) — w^/^ — aw^-^^^ which gives half of an ellipse with axes 1 — a and 1 + a. 

At the special value k = 2, we set k = 2 — S and again take the limit (5 — s- 0, /3 — > with an arbitrary ratio 
a — P/S < 1. The hypergeometric function becomes 1 — a/w^ and the conformal map f{w) = w — a/w gives an ellipse 
with axes 1 — a and 1 + a. 
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